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Abstract — RSS-based device-free localization (DFL) monitors 
changes in tlie received signal strength (RSS) measured by a 
network of static wireless nodes to locate people without requiring 
them to carry or wear any electronic device. Current models 
assume that the spatial impact area, i.e., the area in which a 
person affects a link's RSS, has constant size. This paper shows 
that the spatial impact area varies considerably for each link. 
Data from extensive experiments are used to derive a multi-scale 
spatial weight model that is a function of the fade level, i.e., the 
difference between the predicted and measured RSS, and of the 
direction of RSS change. In addition, a measurement model is 
proposed which gives a probability of a person locating inside the 
derived spatial model for each given RSS measurement. A real- 
time radio tomographic imaging system is described which uses 
channel diversity and the presented models. Experiments in an 
open indoor environment, in a typical one-bedroom apartment 
and in a through-wall scenario are conducted to determine the 
accuracy of the system. We demonstrate that the new system is 
capable of localizing and tracking a person with high accuracy 
(< 0.30 m) in all the environments, without the need to change 
the model parameters. 

Index Terms — wireless sensor networks, device-free localiza- 
tion, tracking 



I. Introduction 

IN recent years, wireless sensor networks (WSNs) have 
been used often for indoor localization. In one of the most 
common approaches, localization is carried out by measuring 
the received signal strength (RSS) of the links composing the 
network. RSS-based localization can be divided into active and 
passive. Active localization is the practice of locating a person 
or asset that is carrying an RF tag. Passive localization does 
not require the person or asset to carry any electronic device, 
sensor, or tag. The division between these two categories 
could also be addressed as follows: in the active case, the 
tracked entity is willing to be located and is co-operating 
with the system, whereas in the passive case the tracked 
object is not co-operating with the system, and possibly 
wants to avoid being located. In this case, the changes in 
the propagation patterns of the wireless network caused by 
movements of people can be exploited for localization. These 
networks are referred to as RF sensor networks HI, since the 
radio frequency (RF) itself is used as the sensing modality. 
This paper focuses on device-free localization (DFL) |T|, also 
referred to as passive localization |2l, sensorless sensing |[3], 
or RF tomography lU, Q. 
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In recent years, RSS-based DFL has emerged as an attractive 
technology for passive indoor localization. The technology, 
which can be used in many applications, such as security and 
surveillance and assisted Uving and elder care, has proven 
its accuracy in open environments f6l, fl], in real-world, 
cluttered environments [2J, |8|, |9|, |10| and also in through- 
wall scenarios ifTTl . lfT2l . An RSS-based DFL system brings 
several advantages over other traditional technologies by being 
able to work in obstructed environments, see through smoke, 
darkness, and walls, and by avoiding the privacy concerns 
raised by video cameras. Moreover, a DFL system can be 
composed of low-cost wireless sensors which measure the 
RSS, in contrast to e.g. ultra-wideband radars, which are also 
RF-based and are capable of measuring the amplitudes, phases, 
and time delays of the multipath signals existing in the radio 
channel, but are also prohibitively expensive. 

Many model-based methods have been already proposed 
to localize and track people using the temporal variations 
of RSS. However, these methods make one or both of the 
following two assumptions. First, the movement of a person 
affects the RSS measurements only when the person is very 
near the line connecting two communicating transceivers |4|, 
fS). Second, when a person's presence is exactly on this line 
between transceivers, which we call the link line, the sensors 
will strictly observe an attenuation fS), Q of the RSS. In 
open environments where line-of-sight (LoS) communication 
among the nodes is dominant and in networks where the 
distance between the nodes is small, both assumptions are 
valid. However, for cluttered environments and longer sensor 
distances, the two assumptions do not apply. In rich multipath 
environments, the RSS of a link can both increase or remain 
unchanged as the link line is obstructed. In addition, as the 
signals propagate via multiple paths from the transmitter to the 
receiver, it is plausible that a person located far away from the 
link line affects a subset of multipath components by reflection 
[131 or scattering [14|, inevitably causing a change in the RSS. 
For these reasons, cluttered indoor environments demand more 
advanced models to characterize the spatial impact area in 
which a person's presence affects the RSS. 

In this paper, we demonstrate that the spatial impact area 
where human-induced RSS changes are measured varies con- 
siderably for each link of the RF sensor network. In addition, 
the spatial impact area is found to depend on the direction of 
the RSS change, i.e., even for the same link, increases and 
decreases of the RSS are observed over different spatial areas. 
We also show that, in addition to a person obstructing the link 
line, RSS changes vary dramatically depending on the radio 
frequency channel. Extensive experiments are conducted in 
two significantly different indoor environments during which 
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more than 26 million RSS measurements are collected. The 
data are used to build a multi-scale spatial weight model 
which describes the relationship between the RSS change of 
a link and the location of the person. The model is built upon 
the concept of fade level, a measure of whether a link is 
experiencing destructive or constructive multipath interference. 
In addition, we derive a measurement model which calculates 
the probability of the person being inside the modeled spatial 
impact area. Real-time algorithms are developed which exploit 
the new models and the performance of the system is vali- 
dated in three different indoor environments, all significantly 
different from one another. In each environment, the methods 
proposed in this paper improve considerably the localization 
accuracy. Furthermore, the more challenging for localization 
the environment, the greater the enhancement in accuracy 
provided by the new methods is. 

II. Related Work 

Several works have already shown that human presence and 
motion alters the way radio signals propagate, enabling the 
localization and tracking of people. Different measurement 
modalities, models, algorithms and applications have been 
proposed, having all as objective to locate people with high 
accuracy. Most approaches achieve sub-meter localization ac- 
curacy. However, a comparison of the results obtained by these 
systems is difficult since they differ considerably in nodes 
number, type of indoor environment, hardware, communica- 
tion protocols, size of the monitored area, just to name a few 
parameters. In the following, we present the characteristics of 
some of the already existing DFL systems. 

One approach to DFL is to estimate the changes in the RF 
propagation field of the monitored area and then to image this 
change field. This image can then be used to infer the locations 
of people within the deployed network. Estimating the changes 
in the propagation field is referred to as radio tomographic 
imaging (RTI), a term coined in [4|. Several measurement 
modalities have been proposed for the purpose of RTI. In 
Q, the attenuation of every voxel in the monitored area is 
estimated using the RSS measurements of many links of the 
network. Attenuation-based RTI is capable of achieving high 
accuracy in small and unobstructed deployments, however, in 
cluttered environments the system loses its capability to locate 
people. 

In ifTTI . variance-based RTI (VRTI) is presented, which is 
capable of localizing and tracking people even through walls 
of a building. A drawback of VRTI is that it is not capable of 
localizing stationary people since the measurements are based 
on a windowed variance of RSS. In addition, VRTI is prone 
to intrinsic motion which is variance of the RSS that is not 
caused by people. The problem is confronted in ifTSll . where 
methods to reduce the noise of VRTI images are presented. 
Results indicate that the tracking accuracy can be improved 
significantly by removing the intrinsic noise of the estimated 
images. 

A measurement modality capable of locating both stationary 
and moving people is presented in |16|. Further, the proposed 
system can achieve high accuracy in open, cluttered, and 



even in through-wall environments. The system is based on 
calculating the kernel distance fYT\ of two RSS histograms, a 
long-term histogram representing the RSS measurements when 
the link line is not obstructed, and a short-term histogram that 
is capable of capturing the temporal RSS variations when the 
person is in close proximity to the wireless link. 

The works in 15], lUS], lHH, EqI also estimate the change 
in the RF propagation field. However, these works derive 
the image from a reduced number of links' measurements. 
In [18;|, distributed processing is applied to minimize the 
information transmitted to a central base station, reducing 
the communication overhead and packet size and increasing 
energy efficiency. In IS), |fT9l . ||20l . the use of compressed 
sensing, i.e., reconstructing a signal from a reduced number of 
measurements is investigated. With this method, the estimated 
images have sharp contrast, while the image formation requires 
fewer links and can be distributed among the sensors of the 
network. However, these approaches have been demonstrated 
only in simple indoor environments and their accuracy in 
challenging multipath environments has not been proven yet. 

The use of channel diversity to enhance DFL accuracy is 
presented in |21|. The work ranks the different frequency 
channels used for communication based on two parameters: 
packet delivery ratio and channel fade level. The results 
demonstrate that fade level is a more important factor than 
communication performance when accurate DFL is required. 
Through channel diversity, the localization accuracy is shown 
to improve by an order of magnitude compared to a sys- 
tem communicating on a single channel. The multi-channel 
system in |21| is used in a long-term residential monitoring 
application in |9J. The work demonstrates the applicability of 
RSS-based DFL in real domestic environments and the ability 
of this technology to achieve high localization accuracy over 
extended periods of time. 

A drawback of imaging-based DFL systems is that they first 
estimate the changed RF-propagation field and then the coor- 
dinates of the person. In this two-step process, information can 
be lost and additional measurement noise can be introduced. 
Hence, methods to estimate the person's location directly from 
the RSS measurements are provided in jS), ifTOl . llT2l . In 
|6|, a particle filter is applied to simultaneously estimate the 
location of the sensors and the coordinates of a person moving 
inside the monitored area. In flOl . a fade level skew-Laplace 
signal strength model and a statistical inversion method are 
introduced to estimate the location of people. In |12|, an 
online learning algorithm is used to determine whether a link 
is obstructed by a person or not, and a particle filter is used 
to locate the person. These works rely on sequential Monte 
Carlo methods to estimate the position of the person. For this 
reason, due to the computational complexity of particle filters, 
especially when the number of used particles is high, these 
systems can not work in real-time. 

In the works presented above, the sensors composing the 
network are deployed around the monitored area. In |22|, the 
sensors are mounted above the monitored area, i.e. hanging 
from the ceiling, and a best-cover algorithm is applied to 
estimate the person location when the RSS changes exceed 
a predefined threshold. The work is extended in ||23| . where a 
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clustering algorithm is used to track two people. 

In this paper, we adapt an imaging based solution f4\, 
Q, ID, Gil, d, GSl, im to estimate the changes in the 
RF propagation field. However, these works assume identical 
spatial areas for links' area of impact, whereas our algorithm is 
the first imaging-based DFL approach to set each link's spatial 
impact area differently. Extensive experiments are conducted 
to derive a multi-scale spatial weight model which more 
accurately describes the human induced RSS-changes with 
respect to the spatial location of people. In addition, a new 
measurement model is introduced which gives the probability 
of the person locating inside the modeled area. The multi-scale 
weight and measurement models are built upon the concept 
of fade level. Channel diversity is exploited to enhance the 
accuracy of the system as in |9|, |21|. However, due to the 
more accurate weight and measurement models, the system 
achieves its best accuracy when all the used channels are 
utilized to estimate the RF propagation field. 

III. Radio Tomographic Imaging 

Radio tomographic imaging is the process of estimating the 
changes in the propagation field of a deployed RF sensor 
network. Several RTI methods have been already proposed 
10 , iflTI . ||T6l. In the following, we describe attenuation-based 
RTI Q, in which the change of the propagation field to be 
estimated is shadowing, i.e., the attenuation of radio signals. 

Let us consider a network that is not occupied by people 
for a certain period of time. During this time, we can measure 
the RSS of the L wireless links. The sample mean of link I 
represents the RSS of the link when this is not affected by 
people. We denote this as f;. We can estimate the attenuation 
of Unk I at time k from the change in received signal strength: 



An(fc) = n(fc) - n, 



(1) 



where r;(fc) is the RSS of link I measured at time k. Link 
I experiences shadowing when RSS decreases (An(fc) < 0). 
To simplify the notation used throughout the rest of the paper, 
time k is left out whenever possible. In attenuation-based RTI, 
the measurement vector y is the change in links' shadowing, 
i.e., y = [Ari . . . Ar^]-'". Measurement yi is a synonym for 
Ar;, however, in Section |IV] we redefine yi. 

In RTI, the attenuation of a link is assumed to be a spatial 
integral of the RF propagation field of the monitored area. 
Some voxels of the discretized area will affect a particular 
link's RSS, and some will not. In our discretized model, this 
corresponds to the fact that each link's change in shadowing 
is assumed to be a linear combination of the change in voxel 
attenuations: 
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where Xj is the change in attenuation in voxel j, aij the weight 
of voxel j for link I, bi the measurement noise, and N the 
number of voxels in the image to be estimated. 

The weighting aij represents a spatial model for how voxel 
j's attenuation affects link Ts measurements. In RTI, it is 
assumed that an elliptical model represents the geometrical 



relationship between the links and voxels. This model, where 
the transmitter and receiver are located at the foci, has been 
effectively applied in g), Q, El, CB], 121J. In this model, 
the voxels j that are within the elUpse of link I have their 
weight aij set to a constant, which is inversely proportional to 
the square root of the Euclidean distance of the communicating 
nodes. Otherwise, the weight aij is set to zero. The weighting 
is mathematically formulated as follows: 
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where d is the distance between the transmitter and receiver, 
d*^ and are the distances from the center of voxel j to the 



transmitter and receiver of link /, and A is the ellipse excess 
path length, a parameter that can be used to tune the width of 
the ellipse. 

When all the links of the RF network are considered, the 
changes in the attenuation field of the monitored area can now 
be modeled as: 

y = Ax + b, (4) 

where y and b are L x 1 vectors representing the measured 
RSS difference and noise of the L links, x is the iV x 1 radio 
tomographic image to be estimated, and A is the weight matrix 
of size Lx N, in which each column represents a single voxel 
of the image and each row the weight of each voxel for a 
particular link. The linear model for shadowing loss is based 
on the correlated shadowing models described in Q, ll24l . and 
on the work presented in 17 |. 

Estimating the image vector x from the link measurements 
y is an ill-posed inverse problem, where the same set of 
link measurements can lead to multiple different images, i.e., 
solutions. Therefore, regularization is required. Here, we use 
a regularized least-squares approach, as in HJ, ifTSl : 
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where af^ is a regularization parameter The a priori covari- 
ance matrix C^. is calculated by using an exponential spatial 
decay: 



3^ 



(7) 



where dji is the distance from the center of voxel j to the 
center of voxel i, is the variance of voxel attenuation, and 6c 
is the voxel's correlation distance. The linear transformation 11 
can be calculated beforehand, once the positions of the sensors 
are known, and it has to be computed only once, enabling real- 
time image reconstruction via (|5]|. 

A. Model Errors 

The linear model for shadowing loss is based upon two 
assumptions. First, a person located on the link line causes 
attenuation. Second, the person affects the RSS measurements 
only when located near the link Une (the value of A is typically 
set to be small Q, |j9l, ifTTl '). Let us validate these two 
assumptions with a test in which 30 nodes communicating on 
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Fig. 1. RSS changes on two different cliannels for link ([0.0 5. 8], [7. 6 4.6]) due to liuman motion. In the figures, the points represent the true location of the 
person when the link measures a decrease, yi < —5 dB, or an increase, yi > 5 dB, of the RSS. The measurements on channel 15 are shown in (a), whereas 
the measurements on channel 26 are shown in (b). 



multiple frequency channels are deployed in an unobstructed 
indoor environment around a square perimeter During the test, 
a person walks inside the monitored area along a predefined 
trajectory. 

In Fig.[T[ the RSS changes of the same link on two different 
channels are shown. In the figure, the points represent the true 
position of the person when large RSS changes are measured. 
For clarity, only large decreases, < — 5 dB (blue squares), 
and increases, yi > 5 dB (red circles), are shown. The 
measurements on channel 15, shown in Fig. [T](a), fit well the 
linear model for shadowing loss. The attenuation of the RSS 
signal is measured when the person is located near the link 
line, thus the area where shadowing is observed is accurately 
predicted by the geometrical ellipse model. However, the 
measurements of the same link on channel 26, shown in Fig. 
[T] (b), indicate that on average an increase in signal strength 
is measured when the person locates inside the ellipse. In 
addition, large RSS changes {\yi\ > 5 dB) are measured 
even when the person is located far away from the link line. 
Consequently, both the measurement and the ellipse weighting 
models are inaccurate for this link and channel. 

IV. Measurement-based Modeling 

In this section we describe the conducted measurement 
campaign that is used to derive the new multi-scale spatial 
weight model and the related measurement model. Section 
|IV-B| provides motivation and insight into why more detailed 
spatial weight models are needed. The section is followed by 
the presentation of the multi-scale spatial weight model and 
the new measurement model. 

A. Data Collection 

To derive more accurate weight and measurement models, 
extensive experiments are conducted in two significantly dif- 
ferent environments: an open indoor environment (experiment 



1), in which all the sensors have LoS communication among 
each other, and a single-floor, single-bedroom apartment (ex- 
periment 2), where multipath propagation is common. 

In both experiments, the sensors deployed in the monitored 
area are Texas Instruments CC2531 USB dongles, set to 
transmit at the maximum power, i.e. -1-4.5 dBm [25 1. The 
sensors communicate in TDMA fashion on multiple frequency 
channels. The IEEE 802.15.4 standard ll26l specifies 16 chan- 
nels within the 2.4 GHz ISM band. They are numbered 
from 11 through 26 and are 5 MHz apart, having a 2 MHz 
bandwidth. The carrier frequency (in MHz) of channel c is: 

/c = 2405 + 5 -(c- 11), ce [11,26]. (8) 

The sequence of transmission is defined by the unique 
sensor's ID number In each packet, the sensors include their 
ID and the most recent RSS measurements of the packets 
received from the other sensors of the network. If a packet 
is dropped, the next sensor in the schedule transmits after 
a back-off time, thus increasing the network's tolerance to 
packet drops. At the end of each communication cycle, the 
sensors switch synchronously to the next frequency channel 
found in a list pre-defined by the user. On average, the time 
interval between two consecutive transmissions is 2.9 ms. Due 
to the low latency between transmissions, the human induced 
changes in RSS of the different channels are correlated. The 
two experiments are described in detail in the following 
paragraphs. 

In experiment 1, shown in Fig. [2] (a), 30 sensors are 
deployed on the perimeter of the monitored area (70 m^). 
The sensors are placed on podiums at a height of one meter 
The sensors are programmed to communicate on channels 
11, 17, 22, and 26 (no Wi-Fi interference). During the test, 
a person walks at constant speed along a rectangular path. 
The trajectory is covered multiple times to collect a sufficient 
number of RSS measurements. Markers are placed inside 
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(a) (b) (c) 

Fig. 2. In (a), an image of the open environment, in (b) the apartment, and in (c) the through-wall scenario is shown. 



the monitored area for the test person to follow, while a 
metronome is used to set a pre-defined walking pace. In this 
way, each collected RSS measurement can be associated to 
the true location of the person. 

In experiment 2, shown in Fig. |2] (b), 33 sensors are 
deployed in a single-floor, single-bedroom apartment (58 m^). 
Most of the sensors are attached on the walls of the apartment, 
while a few of them are placed elsewhere, e.g., on the edge 
of a marble counter in the kitchen or on the side of the 
refrigerator The antennas of the sensors are detached from 
the walls by a few centimeters to enhance the localization 
accuracy |9|. In the experiment, the sensors are programmed 
to communicate on channels 15, 20, 25, and 26 in order to 
avoid the interference generated by several coexisting 802.11 
(Wi-Fi) networks, which would increase the floor noise level 
of the RSS measurements |27|. As in experiment 1, a person 
walks along a pre-defined path at constant speed several times. 

In both environments, the tests are repeated three times: two 
of the tests are used to build the new models presented in this 
work, whereas one of the tests is used to validate the improved 
accuracy obtained by using the new models. To derive the 
new models, more than 26 million RSS measurements are 
collected. Each measurement is associated to the true location 
of the person. Since the new models are based on the changes 
in RSS, an initial calibration at the beginning of each test is 
performed while the monitored area is not occupied by people. 

B. Link Measurements 

People found in close proximity of a wireless link cause 
changes in the RSS by shadowing, reflecting, diffracting, or 
scattering a subset of its multipath components lilTl . 
In open environments in which LoS communication among 
the sensors is dominant, the sensors measure on average a 
decrease in RSS when a person is located on the link line. 
The attenuation of the RSS signal is caused by shadowing. 
This phenomenon has been successfully applied in attenuation- 
based DFL li5J, la, m, GSl, ED- 

In cluttered environments in which multipath is common, 
the changes in RSS due to the presence of people become 
more unpredictable. As the link line is obstructed, the RSS 
can increase, decrease, or remain unchanged ITOl . ifTTl . 11211 . 
In addition, a person can impact the RSS even when located 
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Fig. 3. The RSS measurements of a single link on two channels in different 
fade level conditions. 



far away from the line connecting the sensors |21|. Temporal 
variations of the RSS have been successfully applied to track 
moving people ifTTl and also to locate and track both stationary 
and moving people |jTO| in multipath rich environments. 

The statistics of steady-state, narrow-band fading are related 
to the changes in RSS due to the human presence as described 
in ifTOl . In this work, the authors define the relation using 
fade level, a continuum between two extremes: anti-fade and 
deep fade. A link in a deep fade experiences destructive 
multipath interference and will measure on average an increase 
in RSS when obstructed. In addition, links in deep fade are 
also sensitive to movement far from the link line, measuring 
high variance of the RSS. On the contrary, links in anti-fade 
experience constructive multipath interference. On average, 
these links remain unchanged when people are moving far 
from the link line and measure a decrease in RSS when the 
link line is obstructed. Due to these characteristics, it is the 
links in anti-fade that are the most informative for RSS-based 
DFL, since for them the area where the person affects the RSS 
is small, i.e., a narrow ellipse having TX and RX at the foci. 
Furthermore, for anti-fade links the direction of RSS change 
is predictable, thus easier to model. 

To illustrate the effect of fade level on the temporal RSS 
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Fig. 4. In (a), the combined distribution of the fade levels in experiments 1 and 2. RSS measurements of links that have a fade level of _F £ [7.86, 9.42] 
dB as a function of A are shown in (b). 



changes, the measurements of a single link on two different 
channels are plotted in Fig. [3] The link and channels consid- 
ered in Fig. |3] are the same as the ones shown in Fig. [T] The 
monitored area is empty in the beginning of the experiment 
until a person enters it (k = 185) and starts walking along a 
predefined trajectory. The link line is intersected two times (k 
= 196 and k = 213) within the considered period of time. The 
difference in fade level of the two channels is almost 12 dB. 
The channel in anti-fade (channel 15) measures attenuation 
when the link line is crossed. On the contrary, the deep fade 
channel (channel 26) measures an increase in signal strength 
as the link line is crossed. In addition, the measurements of the 
deep fade channel fluctuate much more than the measurements 
of the anti-fade channel when the person is moving off the link 
line. 

The measurements in Figs. [T| and [3] and previous works 
found in the Uterature HO), IE], El] show that the RSS 
measurements are very different for links in different fade 
levels. This motivates us to further explore how the fade level 
of the links relates to the spatial impact area of the person, the 
measured RSS and its direction of change. In this paper, we 
propose more sophisticated models, derived from experimental 
data, which describe more accurately the area where human 
motion induces RSS changes and the temporal variations of 
the RSS. 

C. Fade Level 

In obstructed environments, where non-LoS communication 
is dominant, the signals travel via multiple paths from trans- 
mitter to receiver, a phenomenon called multipath propagation. 
At the receiver, a phasor sum of the radio waves impinging on 
the antenna determine the strength of the received signal. This 
phasor sum may be constructive or destructive, depending on 
the phase of the waves. This results in a RSS that is a function 
of the center frequency and position in space, an effect called 
multipath fading, which causes a significant deviation from the 



theoretical radio propagation model. The difference between 
the radio propagation model and the mean RSS of link I on 
channel c, f^j, is what we call fade level, and can be estimated 
as follows: 

F,^i^f,,i~P{d,c), (9) 

where P{d, c) is a model for the RSS vs. distance and channel. 

In a wireless network, the signal strength measured at the 
receiver can be modeled with the log-distance path loss model 
HH: 

P(d,c) = Po(c)-10rylogio-f, (10) 

"0 

where Pq{c) is what we call transmit power-normalized ref- 
erence loss on channel c and rj is the path loss exponent. The 
transmit power-normalized reference loss is defined as: 

Po(c) = Pt(c)-Po, (11) 

where Pt{c) is the normalized transmit power on channel c 
and Pq the reference loss at a short reference distance do. The 
normalized transmit power Pt{c) is a linear function of the 
used frequency channel, which must be derived from experi- 
mental data. For the wireless sensor used in the experiments 
f25], the lower frequency channels measure lower RSS values 
than the higher ones, because of the differences in antenna 
impedance matching across a wide frequency band |29|. We 
find the Hnear relationship Pt{c) = (0.1452 -0+ 1.7332) dB 
matched the measured transmit power closely. Normalization 
is required to avoid bias in the fade level estimates. 

To derive the path loss exponent rj in ( [T0| , a calibration is 
performed at the beginning of each test. At the end of this 
period fc,/ for every link and channel is calculated. Since the 
distance between the nodes is known, a least-squares linear 
fit can be used to determine the path loss exponent rj. In an 
indoor environment, it is typical that some sensors measure on 



average RSS lower than what predicted by ( 10 1, whereas other 
sensors have a positive fade level. The combined distribution 
of the fade levels in experiments 1 and 2 is shown in Fig. |4|a). 
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Fig. 5. In (a), measured and modeled A" 
function of fade level. 



and A"*" values as a function of the fade level. In (b), the measured and modeled /3 and /S"*" values shown as a 



D. Multi-Scale Spatial Weight Model 

During the measurement campaign, we determine two mea- 
sures while the person is moving inside the monitored area. 
First, the change in RSS for each time instant fc, Arc,i(fc). 
Second, the minimum excess path length for each link at each 
time instant fc, i.e., A/(fc) = d^j +d'^j —di, so that the person's 
location is on the perimeter of the ellipse. The excess path 
length A; measures the distance of the person from the link 
line. 

We divide the RSS measurements of the experiments in 
20 equally likely bins each containing the same amount of 
measurements based on their fade level. Fig. |4|b) shows the 
change in RSS as a function of A for the links of bin F E 
[7.86,9.42] dB. For each bin we determine the 5th-percentile 
of measurements that account for the largest decreases in RSS. 
For this set of measurements, we derive the median excess path 
length, which we denote as A~ . Correspondingly, we calculate 
A+ for the 5th-percentile of measurements that account for the 
largest increases in RSS. The values of A^ and A+ for the 
different fade level bins are shown in Fig. |5] (a). 

The results indicate that links in a deep fade, i.e. with a 
negative fade level, measure decreases in RSS within a large 
area. The excess path length of the ellipse decreases as the fade 
level increases, confirming the fact that, when RSS attenuation 
is observed, anti-fade links are more informative for DFL. 
The opposite is true when an increase in RSS is measured. 
In this case, links in a deep fade measure the human induced 
RSS change within a smaller area than the anti-fade links. 
However, the difference in A+ for deep fade and anti-fade links 
is relatively small compared to the difference in A^ values. 

An exponential decay model accurately fits the A~ values. 
Correspondingly, an exponential growth model is used to fit 
the A+ values. For both, a least-squares fit is used to derive 
the model parameters. Now, A can be determined separately 
for each link, frequency channel and direction of RSS change 



TABLE I 
Ellipse model parameters 



Parameter 


<5 is — 


5 is + 




-5.7874 


102.7284 


K 


0.2112 


0.5016 



based on the fade level of the link: 

(12) 

where 5 indicates the direction of RSS change, i.e., (5 is — for 
measured decreases in signal strength and + for the measured 
increases. Table [l] shows the values of the parameters of the 
ellipse model. 

E. Measurement Model 

In the previous section it was shown that the area in 
which a person affects the link's RSS measurements strongly 
depends on the fade level of the link and the direction of 
the RSS change. In this section, we derive a measurement 
model to determine the probability of the person being within 
the area predicted by the new multi-scale weight model. The 
new measurement model is based on both the magnitude and 
direction of RSS change, and fade level of the link. 

The RSS measurements are organized as in the previous 
section. This time, the measurements of each fade level bin 
are further divided into 100 bins so that each bin contains 
Arc,; measurements having approximately the same value. 
For each bin, we estimate the probability of the person being 
outside the ellipse. Figure |6|a) shows the probabilities of the 
bin Fi e [7.86, 9.42] dB (i.e., anti-fade links) that the person 
is located outside the modeled ellipses for each Ar^;. The 
negative x-axis of the figure shows the probability of the 
person being outside the ellipse A^ when a decrease in signal 
strength is measured. Correspondingly, the positive x-axis 
shows the probability of the person being outside the ellipse 
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A+ when an increase in RSS is measured. The probabilities 
of bin Fi e [-13.70,-10.13] dB (i.e., deep fade links) are 
illustrated in Fig. |6|b) for comparison. 

The probability that the person is located inside the modeled 
ellipse depends on the magnitude of Ard - The larger the RSS 
change is, the more likely it is that the person is located inside 
the modeled ellipse. It is again the anti-fade links that are 
more informative and trustworthy. For example, when a link 
measures a 10 dBm attenuation, the anti-fade link (Fi — 8 dB) 
has a 69% probability that the person is inside the modeled 
ellipse. In comparison, the probability is the same with the 
deep fade link (Fi — —8 dB). However, for the deep fade link 
the modeled ellipse is considerably larger, i.e., A+ — 0.8413 
whereas = 0.0530. Thus, we can also state that the anti- 
fade links indicate the person's location more accurately when 
increases in signal strength are measured. For example, when 
a link measures a 10 dBm increase in RSS, the anti-fade link 
(Fi — 8 dB) has a 97% probability that the person is inside 
the modeled ellipse. In comparison, there is a 63% probability 
with the deep fade hnk (F/ = -8 dB). 

Different models are used for the negative and positive RSS 
changes. The decay rates /3~ for the different fade levels when 
Arc,; < are nearly constant thus a fixed value is used in 
( 14 1. On the contrary, decay rate /3+ is a function of the fade 
level. Thus, /3^; can be determined separately for each link 
and frequency channel based on the fade level of the link as 
follows: 



TABLE II 

Measurement model parameters 



(13) 



The measured and modeled decay rates are shown in Fig. 
[5jb) and the parameters of the models are listed in Table [ll] 
Now, the probability of the person being located inside the 
modeled ellipse at time k and measured RSS change Ar^^i 
can be estimated as: 



(14) 



/9- 






0.1172 


13.0018 


0.1839 



V. Methods 

The image reconstruction procedure for RTI can be used 
as a theoretical framework for estimating the changes in the 
RF propagation field with the new weight and measurement 
models. However, minor adjustments need to be made to 
RTI as it was introduced in Section [III] First, instead of 
applying the changes in RSS as given in ([T]), we apply the 
probability of the person being located inside the modeled 



ellipse which is given in ( 14 1. Thus, the new measurement 
vector on channel c when attenuation is measured is y~ = 
[p^]^, . . . ,p^^]. Correspondingly, for measured increases in 
RSS y+ = [p'^n ■ ■ ■ tP^l]- Thus, the complete measurement 
vector on channel c becomes yc = [ycMycT]- When consid- 
ering all the channels, the complete measurement vector is 
y — [yi| • ■ ■ lyc]"'", where C is the number of used channels. 

The spatial weighting model in Q has to be reformulated 
since A is now unique for each link. The new model can be 
mathematically expressed as follows: 




if d\^+di- <d+\ii 

otherwise 



(15) 



where ; ^ is the weight of voxel j for link I on channel c 
for RSS change direction 5. Because the area covered by the 
ellipses varies, we weight less the links that cover a larger area 
by setting the weight to be inversely proportional to the area 
of the ellipse, i.e. Uj -p^, where rij is the number of voxels j 
that are within the ellipse of link I and p^ is the area of the 
voxel. 

Now, when all the links and channels of the RF network 
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are considered, the changes in the RP propagation field of the 
monitored area can be modeled as: 



TABLE III 

Image reconstruction parameters used in the experiments 



y = Wx + n, 



(16) 



where y and n are the measurement and noise vectors, 
correspondingly, both of size 2LC x 1. As in Q, x is the image 
to be estimated, and W is the new multi-scale weight matrix, 
of size 2LC x N. The regularized least-squares approach in 
(j6]l can be used also with the new models by substituting the 
weight matrix A with the new weight matrix W. 



A. Localization and Tracking 

From the estimated image (|5]l, the position of a person can 
be determined by finding the voxel of the image that has the 
maximum value: 



J = arg max x. 

N 



(17) 



The position estimate is therefore z — Zj, where Zj represents 
the coordinate of the center of voxel j. For tracking the 
movements of the person, we apply a Kalman filter ll30l 
designed to estimate the position, velocity, and acceleration 
of the tracked person as in [18]. Note that the localization 



proposed in (17 1 is only capable of locating one person. Multi- 
target localization and tracking is a more challenging issue in 
DFL, and it is outside the scope of this paper Thus, we do 
not propose coordinate estimators for the multi-target case. 
Readers are referred to flOl, 123], IHl, lEH for RSS-based 
multi-target tracking. 

In Section |VT] we evaluate the accuracy of the system. For 
that purpose the localization error at time k is defined as the 
distance between the estimated and true position: 

efe = \\zk - ZH,kh (18) 
where zn^k is the true human location. 

B. Performance Benchmarking 

The system presented in this paper, which we denote as 
multi-scale RTI (msRTI), is compared to two other systems: 
an RTI system that ranks the channels based on the fade level 
II2TI . ||9l and a RTI (see Section IIIi system that exploits 
channel diversity. The fade level RTI (flRTI) system ranks the 
different channels based on their fade level and uses only the 
measurements of the m most anti-fade channels. The RSS 
measurements collected on the selected channels are averaged 
to obtain the measurement vector y. The channel diversity RTI 
(cdRTI) system considers each channel as an independent link, 
i.e., as if at each node's location there were multiple sensors 
communicating on different channels. For the two systems 
used for benchmarking, the image reconstruction parameters 
that maximize the localization accuracy are derived through 
simulations. For the msRTI system proposed in this paper, the 



Parameter 


Value 


Description 


P 


0.1524 


Voxel width [m] 


O'x 


0.0316 


Voxels standard deviation [dB] 




1.4142 


Noise standard deviation [dB] 


Sc 


4 


Decorrelation distance [m] 



image reconstruction parameters, shown in Table III are the 
same for all the tests. 



C. Experiment description 

In the beginning of the tests (see Section IV-Ai, a one 
minute calibration period is performed. During this period, 
the mean RSS and fade level of every link and channel are 
calculated. Online training and calibration of the network are 
not within the scope of this paper. However, both are research 
topics that will be addressed in future work. For now, readers 
are referred to ifTTl . lfT2l . ||9l for systems that do not require 
calibration {TT\, have online training lfT2ll . or are capable of 
adapting to the changing environment ||9l- 

Experiment 1 is solely a tracking experiment. In the test a 
person enters the monitored area, after the initial calibration, 
and walks along a predefined rectangular path, covering 14 
laps before leaving the area. In experiment 2, both localiza- 
tion and tracking accuracy are evaluated. In the localization 
experiment, a person is standing still at one of 15 known 
positions for a predetermined period of time so to compute 
multiple position estimates at each location. The 15 positions 
are chosen so as to cover all the areas of the apartment. During 
the experiment, each position is visited multiple times. In the 
tracking experiment, a person is moving at constant speed 
along the paths connecting the same 15 locations used in the 
localization experiment. 

Since the models presented in this paper are derived from 
the data collected in the environments of experiments 1 and 2, 
it is important to validate their generality in an environment 
differing from the ones used to build them. Hence, we conduct 
an additional test in a more challenging through-wall scenario 
(70m^). In the experiment, 30 nodes are deployed around a 
lounge room, outside the walls surrounding the monitored 
area. The sensors are set on podiums as in experiment 1 . Five 
channels (11, 15, 18, 21, and 26) are used for communication. 
During the experiment, the person is standing still at one of ten 
predefined positions, located along a rectangular perimeter so 
to evenly cover the entire monitored area. Besides evaluating 
the localization accuracy of the system in a through-wall 
scenario, we also focus on analyzing the effect that the number 
of frequency channels used for communication has on the 
accuracy. 

VI. Experimental Results 

A. Experiment 1, open environment 

The true and estimated trajectory of the person during the 
test are shown in Fig. |7|a). The coordinate's estimates of two 
consecutive laps of the test are shown in Fig.|7|b). The msRTI 
system is capable of tracking the person with an average error 
of e = 0.17 m, whereas the average error for flRTI and cdRTI 
systems are e = 0.22 m and e — 0.39 m respectively. In 
ETI . we demonstrate that an RTI system exploiting channel 
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Fig. 7. In (a), the estimated and true trajectories of the person in the open environment test and in (b), the coordinate estimates of the person over two 
consecutive laps. In (c), CDFs of the tracldng accuracy with the tliree different methods in the open environment experiment. 
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Fig. 8. In (a), the estimated and true locations of a stationary person. In (b), the estimated and true trajectories of a moving person. In (c), CDFs of locating 
a stationary and moving person with msRTI in the apartment experiment. 



diversity always outperforms an RTI system relying on single 
channel communication and that the performance can be 
increased even further when the channels are ranked based on 
their fade level. However, the drawback of the method used 
in the flRTI system is that the RSS measurements collected 
on the deep fade channels are actually not considered at all. 
In fact, when the traditional spatial weight model is applied, 
not considering these noisy data improves the localization 
accuracy. Instead, when the new methods are applied, the 
RSS measured on all the channels can be effectively used to 
further improve the localization accuracy, as demonstrated by 
the results. 

The cumulative distribution functions (CDFs) of the local- 
ization errors for the three different systems are shown in 
Fig.|7|c). In this experiment, the median error is 0.16 m with 
msRTI, 0.20 m with flRTI and 0.30 m with cdRTI. One major 
advantage of using msRTI is that all the position estimates are 
accurate — the most inaccurate being 0.63 m. With the other 
two systems, some position estimates are more than a meter 
off the true position of the person. The highest errors for the 
two methods are 2.03 m with flRTI and 2.12 m with cdRTI. 

B. Experiment 2, apartment 

The true and estimated positions of the localization ex- 
periment are shown in Fig. [8ja). Respectively, the true and 
estimated trajectories are shown in Fig. |8jb). The average 
localization accuracy of the experiment is e = 0.26 m and the 



average tracking accuracy of the experiment is e = 0.23 m 
with msRTI. The CDFs of the two experiments are shown in 
Fig. [Sjc). The median error of the localization and tracking 
experiments are 0.15 m and 0.18 m, respectively. Corre- 
spondingly, the 95th-percentile of the two errors are 1.24 
m (localization) and 0.59 m (tracking). The location found 
in the bedroom degrades considerably the accuracy of both 
experiments. However, since in the localization experiment the 
time spent in that particular location is longer, the localization 
results are affected more than the tracking ones. 

The combined error of the localization and tracking experi- 
ments is 0.25 m with msRTI, whereas it is 0.26 m with flRTI 
and 0.36 m with cdRTI. The difference in accuracy among the 
systems is not as remarkable as in the open environment. In the 
apartment, the sensors are deployed in all the rooms, and there 
exist many sensors that have LoS communication with each 
other In addition, the distance separating some of the sensors 
is small. Due to this, the links measure on average attenuation 
as the link line is obstructed. Therefore, an attenuation-based 
system that exploits a fixed A model is capable of performing 
adequately well in the apartment deployment. 

C. Experiment 3, through-wall 

The true and estimated positions of the experiment are 
shown in Fig.|9|a). The figure shows for comparison the results 
obtained with flRTI. The average localization error is e = 0.30 
m with msRTI, whereas with flRTI the average localization 
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Fig. 9. The estimated and true positions of the person with msRTI and flRTI is shown in (a). The effect of channel number on the localization accuracy of 
the systems presented in (b). 
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Fig. 10. The estimated change in the attenuation field with four different methods. In (a), using RTI on a single channel. In (b), using cdRTI (each channel 
is considered as a unique link). In (c), the channels are ranked based on their fade level and the three most anti-fade channels are used (flRTI). In (d), with 
the new models presented in this paper (msRTI). The true location of the person is represented by the white cross. 



accuracy is e = 0.70 m, and e = 1.57 m with cdRTI. For 
clarity, the resuhs of cdRTI are omitted from Fig |9]^a). The 
improved performance of the new models is due to the fact that 
in the through-wall scenario the number of links measuring 
attenuation when the link line is obstructed is considerably 
smaller than in the other two environments. Furthermore, in 
this type of environment it is common that human-induced 
changes in the RSS are observed also elsewhere than on the 
link line. This further motivates the use of the multi-scale 
spatial weight model in cluttered environments. 

As demonstrated in ETll . channel diversity enhances the ac- 
curacy of DFL considerably, since it increases the probability 
that at least one (or more) of the channels will be in anti-fade, 
thus favoring the localization attempt. Figure |9|b) shows the 
average localization errors over the ten predetermined loca- 
tions obtained with msRTI for all the possible combinations 
of the five channels used in the through-wall experiment. The 
average localization error obtained with msRTI is smaller than 
the ones obtained with flRTI. The msRTI system achieves a 
better average accuracy by using only two channels than the 
fIRTI by using all five channels. 

The advantage of exploiting channel diversity and the new 



models is shown in Fig. 10 The figure shows the estimated 
change in the propagation field with four different methods. In 



(a), the change in the attenuation field is estimated by using a 
single channel for communication (RTI). With this method, the 
localization error is 5 m. In (b), channel diversity is exploited 
(cdRTI). Even though the position estimate obtained with this 
method is accurate, the image is very noisy, containing two 
separate "blobs" which could be interpreted as two different 
individuals. In (c), the change in the attenuation field is 
estimated by ranking the channels based on their fade level 
and using only the three most anti-fade channels (flRTI). The 
system locates the person with a 0.34 cm error. Not only is 
the image more accurate, it also contains much less noise then 
in (a) and (b). In (d), the new models are applied (msRTI). 
The true and estimated positions overlap, and the image is not 
corrupted by noise. 

VII. Conclusion 

In this paper, we present new methods to enhance the 
accuracy of RSS-based DFL. The improvements concern four 
aspects: exploiting channel diversity, deriving a more accurate 
spatial model for the human-induced RSS changes, proposing 
a measurement model that determines the probability of the 
person being inside the modeled area, and taking into con- 
sideration the direction of the RSS change. In the paper, we 
experimentally show that a person located on the link line not 
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only causes attenuation of the signal but can also cause the 
signal strength to increase. In addition, we show that a person 
can influence the link measurements even far away from the 
link line. Both results are found to depend on the fade level 
of the link— a measure indicating if the link is experiencing 
destructive or constructive multipath interference. The fade 
level is then used to derive more accurate measurement and 
spatial models for the human-induced changes in the RSS. 

The performance of the new models proposed in this paper 
are validated in three different indoor environments (open, 
cluttered, through- wall). The results demonstrate that the new 
method outperforms the current state-of-the-art methods in [7], 
Ell . Moreover, the improvement in localization accuracy with 
the new system is more consistent the more the environment 
is challenging for RSS-based DFL. The results indicate that 
the system here presented is the most accurate RSS-based 
DFL system reported to date, capable of achieving 0.30 m 
localization accuracy even in through- wall scenarios. 

In future work, we will investigate spatial propagation 
models that can adaptively learn the geometrical propagation 
patterns of the multipath components and their associated 
amplitudes. Adapting the models to the environment can 
provide the means to further improve the accuracy of DFL. 
So far, only the spatial correlations of the link measurements 
have been considered when estimating the changes in the 
propagation field. However, the link measurements are also 
time correlated. Future research will explore how to exploit 
the time correlation of the measurements in the context of 
DFL. 
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